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Abstract 



^^ ' In this work we solve the nonhnear second order differential equa- 

^^ tion of the simple pendulum with a general initial angular displacement 

{9{0) = Oq) and velocity (^(0) = ^o), obtaining a closed-form solution 

,S^ , in terms of the Jacobi elliptic function sn(u, k), and of the the incom- 

Qh' plete elliptical integral of the first kind F{(p,k). Such a problem can 

^ I be used to introduce concepts like elliptical integrals and functions to 

C/2 , advanced undergraduate students, to motivate the use of Computer 

^ ' Algebra Systems to analyze the solutions obtained, and may serve as 

Z^ I an exercise to show how to carry out a simple generalization, taking 

jyj . as a starting point the paper of Belendez et al [8], where they have 

O I considered the standard case 9{0) = 0. 

C/3 . Keywords: Simple Pendulum, Large-scale Period, Elliptic Integrals. 

, ^' 1 Introduction 



An usual way to present the simple pendulum model in a first course of 
general physics is to consider a small angular displacement, and in this way 
fSj I to linearize the corresponding second order differential equation [2] . This first 

O ■ approximation lead us to conclude that, for small amplitudes, the period of 

. . oscillation is independent from the amplitude: 

o 
o 




(1) 



In order to present a more realistic expression for the period in the case of 
larger angles, it was proposed another approximated formula, derived from 
the complete elliptic integral of first kind [3l d] : 



T « 27r ./ — ^ (2) 
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where ^o is the initial angular displacement and the amplitude of oscillation, 
supposing the the pendulum starts at rest. 

In fact, series approximations are common in classical mechanics text- 
books [5], and another approximations were recently proposed [6, 7J. 

It is well known that the exact formula for the large-angle period of a 
simple pendulum, involving the elliptic integral of the first kind i^(-lj is 
given by p[]: 

n / a^\ 

(3) 

Belendez et al [8j, in a very didactic note, have derived this formula from 
the integral expression of 6{t). 

In this brief work we will essentially solve the same problem treated 
by Belendez et al [8], but considering a non-zero initial angular velocity. 
We present a closed- form expression for the angular displacement 0[t), in 
terms of the Jacobi elliptic function sn(it, k), and of the incomplete elliptical 
integral of the first kind F{Lp, k). We also made some energy considerations, 
and derive an expression for the period of oscillation T{9o,(I)q), involving 
K{-). We believed that such an exercise can be fruitfully used in order to 
introduce such concepts as elliptical integrals and functions, to encourage 
the use of Computer Algebra Systems (Maple®, Mathematica® , etc.), and 
to show how a simple generalization can be carried out. 

2 The Model 



The dynamics of an ideal simple pendulum is given by the following initial 
value problem: 

0{O) = 9o (4) 

I ^(0) = 00 



+ LOn sin 9 = 



where cjq = w| is its natural angular frequency, 9o £ [0,27r] is the initial 
angular displacement, and (f)Q £ ?R. is the initial angular velocity. 

Following |8j, we start to solve this second order differential equation, 
multiplying it by -^i 

d9£9 (19 o . ^ ^ 

dt dt^ dt " 



and noting that it can be rewritten as: 



di 



2\dt) 



loq cos ( 



^K{-) is implemented in most computer algebra systems available in the market. 



Then, integrating this equation on the time interval [0,t], t > 0, we obtain: 



1 (<^^\ 2 n 



1 fd9 



2\dt 



1 UO 



2 \dt 



t=o 



ujQ{cos9{t) — cos^o 



f — j = 2u;^(cos e{t) - cos 6*0) + 



Now, the objective is to write this equation in order that, when we integrate 
it, we win obtain an ehiptic integral. To do this, we define the new variables: 



y{t) = sin 



e{t) 



1 ■ 2 I ^0 
k = sm — 
2 



and rewrite the equation as: 

2 



(I) ^.S.-(:-^)u-.^).^U-.^ 



where we have used that 

dy dy dO 1 



dt dO dt 



■ cos 



e\ d9 



2 dt 



and that 



Defining 



follows that: 



cos 



1 - sin^ 



r = LdQt 

= JL 



dz dz dy dt 1 dy 

dr dy dt dr loq\'% dt 



Substituting -^ on ([S]) yields: 



dz 
d^ 



[l-,^)^l-kz') + -P^{l-kz^) 



Au'^k' 



(5) 



(6) 



(7) 



= {I - kz'){ji - z') 
where 79 = 1 + t^; ^ ^ (0) 1)) ^^^d the transformed initial conditions are: 

z{0) = 1 



dz 




(t>o 


(Of, 
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cot 




dT 


r=0 


2a;o 


V^ 



Taking the square root on both sides of ([7]) , and after some algebra, we have 
that 

dr = ± , ^" 

Integrating on the interval [0, r] we obtain: 

, 1 /"^ dz 



Vk 



' ^{i-^')hi-^') 



which is equivalent to: 



Vk 



dz 



' V(i--')(^o-^') ^° ^J i-k - ^') i^i - ^') 



dz 



Now, we can write these two integrals in terms of the inverse Jacobi elliptic 
function and of the incomplete elliptical integral of the first kind (Formula 
219.00, pp. 58, 0): 



lb r = sn 



<^ sn 



-■(^.*)-- 



arcsin I — 1 ,7o^ 

70 



(8) 



70 



->7o^ 



F I arcsin I — I , 7o ^ 
V V7o/ 



±T 



Applying the Jacobi elliptic function sn(-, •) on both sides of this equation 

— = sn ( F ( arcsin ( — ) , 7o fc I ± r, 79 A; ) 
70 V V V70/ / / 

and coming back to the original variables 

= 7oVk sn I F I arcsin I — J , j^k J it r, 79 A; 



sm 



we finally obtain the solution for the angular displacement as a function of 
time: 



6{t) = 2 arcsin 



70 



vk sn I F I arcsin I — ] , 79^; 1 ± ujQt, 7q A: ] 



where 



7o = 1 + 



Au'^k 



, k = sin^ 



and LOn 



(9) 



(10) 



Observe that expression ([9]) represents two possible solutions: there are a 
plus and a minus signal in front of ojQt. Then, we will choose the "positive 



branch" if the initial angular velocity is positive {(po > 0), and the "negative 
branch" if the initial angular velocity is non-positive {(po < 0). 
In particular, (J)q = => 7q = 1, and ([9]) reduces to: 

6{t) = 2 arcsin vk sn ( F ( — , fc j — uot, k 

From Formula 110.06 (pp. 9, [9J), the following equality holds: 

F[^,k)=K{k) 

where K{-) is the complete elliptical integral of the first kind. In this way 
we recover the solution obtained by ^ for the case of a null initial angular 
velocity: 

9 (t) =2 aicsin Vksn{K{k) -uJot,k) (11) 

In Figure 1 we plot the solution ([m— lll|) for 6q = 7r/2, using the software 
Mathematica® , and in Figure 2 the corresponding angular velocity cl){t) = 
9{t). The parameters used were: ^o = f rd, / = 1 m, and g = 9.8 m/s^ (these 
same values for / and g will be used in Figures 3 and 4 below). 

3 Some Energy Considerations 

The problem that we are analysing is conservative, because the pendulum 
is only under the influence of Earth's gravitational field, that is, we are not 
considering any kind of dissipative forces. 

It is straightforward to show that the kinetic energy in the system at 
time t is given by: 

Kit) = \mfcp\t) 

where (pit) = 6{t). The potential energy, considering the lowest level of the 
pendulum as the reference point, is: 

U{t) =mgl{l-cos{e{t))) 

Then, the sum of these two quantities is constant over time, and equals the 
mechanical energy, E^ of the system: 

E = mgl{l - cos{e{t))) + -ml^(t?{t) 

In particular, we know the initial angular displacement (^o) and velocity 
{(po). Consequently, the mechanical energy is given by: 



E = mgl{l - cos(6'o)) + -ml'^cpl 



and the following expression must be true for all time: 

5/(1 - cos(0(t))) + h^4>\t) = gl{l - cos(0o)) + ^1^4- 

Dividing this equation by P and multiplying by 2, we can rewrite it as 

2u;^(l - cos{9{t))) + (P^t) = 2wg(l - cos(0o)) + 4>l{t) (12) 



From (|12p . the angular velocity of the pendulum for the some displace- 
ment 9{t) is given by: 

Ht) = ±^2Lo^{cos{9{t))-cos{9o))+(Pl (13) 

3.1 Stopping at the top 

For instance, if we want to know what initial angular velocity 4>q is necessary, 
given that we know 9q, to the pendulum stops at its maximum height, we 
need only to consider 9{t) = tt and cj}{t) = in p^ . obtaining: 



(Ao = (pc := ±a;ov/2(l + cos(eo)) (14) 

In conclusion: if \(Pq\ < \(pc\, the pendulum will oscillate forever around its 
equilibrium point (its lowest level); if |0o| = \4'c\, the pendulum stops at its 
maximum height; and if |(/'o| > \4'c\ the pendulum will execute a circular 
motion in relation to its fixing point. In Figure 3 we show a graphic with 
the relation \(pc\ versus ^o- 

3.2 Pendulum's Period for |0o| < \(pc\ 

For (po = 0, remember that the pendulum's period is given by ([3]): 

where Tq = 27r^/- is the period of the pendulum for small oscillations ([T]). 
If ^0 / 0, we can still use this formula, just replacing 9q by 9max, that is, 
the maximum angular position where the pendulum stops. We can obtain 
this value setting (j){t) = in ([13]), and isolating 9 = 9max'- 

9max = arccos { cos(6'o) - -4>l \ (16) 

Then, the period of the pendulum, for any initial position and velocity (con- 
sidering \(J)q\ < |</>c|), is given by: 

I o /arccos fcosf^o) — ^0n, • • , ^ 

T(^o, M = -ToK W ^ \ " '^°^ (17) 





In Figure 4 we plot T(0O)0o) and verify that the pendulum's period is an 
increasing functions of both 6*0 and 4>o . 

6 



4 Concluding Remarks 

In this exercise we have derived a closed-form solution for the angular dis- 
placement of a simple pendulum in terms of the Jacobi elliptic function 
sn(n, k), and of the incomplete elliptical integral of the first kind F{ip, k): 



9{t) = 2arcsin 
where 



70 



yk sn I F I arcsin ( — j , 79^ ) ± uj^t, 79 /c ) 



2 _ 1 , 't'l t, _ .^„2 / f^o \ _ , ,2 _ 19 



70 = 1 + ^, A: = sin ^- j , and a;o ^^ 

and ^0; 00 ^-re the initial angular displacement and velocity, respectively. 

In addition, we have shown that, for an initial angular velocity satisfying 
|(/>o| < U!Q\J2{1 + cos(0o))) the pendulum will oscillate with a period equals: 

2 / r, / arccos (cos(0o) — ^0n) \ \ 
T(^o, 0o) = -ToK f sin2 f ^ ^ "^ ^^^' j j 

Finally, such an exercise may be used to initiate advanced undergraduate 
students to concepts such as elliptical integrals and functions, to the use of 
Computer Algebra Systems (Maple®, Mathematica®, etc.), and to show 
how a simple generalization can be carried out. 
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Figure 1: Angular Position versus Time 
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Figure 2: Angular Velocity versus Time 
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Figure 3: \(j)c\ versus 9q 




Figure 4: r(6'o,' 
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